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ABSTRACT 


Scalar  electrical  conductivities  have  been  computed  for  nitrogen  and 
argon  plasma  in  thermal  equilibrium.  Calculations  have  been  made  for 
these  gases  both  unseeded  and  seeded  with  cesium,  potassium,  and  NaK 
in  amounts  ranging  from  0.  01  to  10  percent  seed  by  weight,  and  cover  a 
range  of  pressure  and  temperature  from  10-2  to  10^  atm  and  from  2000  to 
15,  000°K,  respectively.  Calculated  electrical  conductivities  are  pre¬ 
sented  in  figure  form  as  a  function  of  temperature  with  seed  rate  and  pres¬ 
sure  as  parameters.  Some  effects  of  applied  magnetic  fields  up  to  20,000 
gauss  on  the  conductivity  of  potassium-seeded  nitrogen  are  also  given. 
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NOMENCLATURE 

Magnetic  flux  density,  weber/m^ 

Ion-slip  factor  defined  following  Eq.  (24) 

Seeding  rate  factor  defined  by  Eq.  (13) 

Electric  field,  volt/m 

Effective  electric  field  (E  +  V  x  B),  volt/m 
Electronic  charge,  1.602  x  10~19,  coulomb 
Fraction  of  neutral  particles,  (1  -  <£) 

Planck's  constant,  6.62  x  10“^,  joule-sec 
Ionization  potential,  joule 
Unit  vectors 

Current  density,  amp/ m2 
Equilibrium  constant,  l/m^ 

Boltzmann  constant,  1.380  x  10“23^  joule/°K 
Molecular  weight,  kg/(kg-mole) 

Particle  mass,  kg 

Avogadro's  number,  6.025  x  10^6 

Number  density,  number/m^ 

Pressure,  newton/m^  or  atm  (1  atm  =  1.013  x  10^ 
newton/ m2) 

Collision  cross  section,  m^ 

Gas  constant,  8.314  x  10^,  joule/ (kg-mole)  °K 

Seeding  rate  (ratio  of  weight  of  seed  to  total 
mixture  weight) 

Temperature,  °K 

Velocity,  m/  sec 

Mean  particle  velocity,  m/sec 

Molar  fraction 

Mass  fraction 

Partition  function 

Defined  by  Eq.  (22) 


IX 


AEDC-TDR-64-1  1  9 


€0  Permittivity  of  free  space,  lCT9/367r,  farad/m 

k  Defined  by  Eq.  (21) 

A  Ratio  of  Debye  length  to  impact  parameter  defined 

by  Eq.  (3) 

a  Scalar  equilibrium  electrical  conductivity,  mho/ m 

<7eff  Effective  tensor  conductivity,  mho/m 

r  Mean  free  time  between  collisions,  sec 

r  Effective  mean  free  time  defined  after  Eq.  (24),  sec 

0  Fraction  of  ionized  particles 

u  Cyclotron  frequency,  1/sec 


SUBSCRIPTS 

A 

Argon  or  applied 

Cs 

Cesium 

e 

Electron 

eff 

Effective 

ei 

Electron-ion 

ej 

Electron- specie 

en 

Electron-neutral 

f 

Fully  ionized 

g 

Gas 

H 

Hall 

i 

Ion  or  induced 

in 

Ion-neutral 

j 

Specie  index 

m 

Mixture  after  ionization 

N 

Nitrogen  atom 

n2 

Nitrogen  molecule 

n 

Neutral 

s 

Seed 

T 

Total 

w 

Weakly  Ionized 

x,y,  z 

Orthogonal  co-ordinates 
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1.0  INTRODUCTION 

Because  of  recent  research  into  magnetohydrodynamic  (MHD)  power 
generation  and  plasma  acceleration,  there  has  been  an  increasing  effort 
to  determine  the  electrical  conductivity  of  gaseous  plasmas.  A  knowledge 
of  the  conductivity  is  essential  for  the  interpretation  of  MHD  experimental 
phenomena  or  for  making  realistic  predictions  from  theory.  Calculated 
values  of  scalar  d-c  electrical  conductivities  are  available  for  a  few  gases 
(see  Refs.  1  through  6  for  example).  Also,  several  investigators  have 
attempted  to  measure  electrical  conductivities  but  usually  over  ©.  limited 
temperature  range  and  below  2500°K.  References  7  through  15  summa¬ 
rize  most  of  the  available  experimental  data  which,  in  general,  have  been 
obtained  for  air  and  seeded  combustion  gases. 

As  part  of  an  experimental  investigation  of  a  steady-flow,  crossed- 
field,  d-c  accelerator  at  Arnold  Engineering  Development  Center  (AEDC), 

Air  Force  Systems  Command  (AFSC),  additional  computations  of  equilib¬ 
rium  scalar  conductivities  have  been  made  for  nitrogen  and  argon  both  un¬ 
seeded  and  seeded  with  cesium  (Cs),  potassium  (K),  and  NaK  (a  eutectic 
mixture  of  sodium  and  potassium).  The  results  of  these  calculations  which 
cover  a  temperature  range  from  2000  to  15,  000°K  and  pressures  from 
0.01  to  10  atm  are  presented  here.  For  unseeded  argon,  however,  the 
results  of  some  additional  calculations  are  given  for  pressures  up  to  10^  atm. 

Application  of  a  magnetic  field  to  a  plasma  gives  rise  to  electric  fields 
and  currents  in  various  directions;  for  example,  the  Hall  field,  and  as  a 
result  the  conductivity  becomes  a  tensor  quantity.  Sonnerup  analyzes  this 
effect  (Ref.  16)  and  Brunner  (Ref.  1)  has  made  some  calculations  to  show 
the  significance  of  tensor  conductivity  in  air.  Some  additional  calculations 
of  the  three  primary  tensor  components  of  the  conductivity  induced  by 
magnetic  fields  from  5000  to  20,  000  gauss  (that  is,  0.5  to  2.0  weber/m2) 
in  seeded  nitrogen  are  also  presented. 

Application  of  an  electric  field  to  a  plasma  can  result  in  the  electrons 
being  selectively  heated,  and  the  conductivity  becomes  a  function  of  the 
applied  field.  Kerrebrock  (Refs.  17  and  18)  has  termed  this  phenomenon 
"non- equilibrium  conductivity.  "  No  attempt  has  been  made  to  account  for 
this  effect  in  any  of  the  calculations  reported  here. 

2.0  SCALAR  ELECTRICAL  CONDUCTIVITY 
2.1  CONDUCTIVITY  EQUATIONS 

Theories  for  predicting  electrical  conductivities  in  gases  have  been 
developed  for  the  two  extreme  cases  of  weakly  ionized  gases  and  fully 

Manuscript  received  May  1964. 
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ionized  gases.  In  both  cases  the  analysis  assumes  thermal  equilibrium 
(Tn  =  Ti  =  Te)  and  single  ionizations  (ne  =  ni). 


Chapman  and  Cowling  (Ref.  19)  have  made  a  careful  analysis  for 
weakly  ionized  gases  in  which  they  consider  binary  collisions  between 
neutral  particles  and  electrons;  that  is,  since  ni  <<nn  they  only  take 
close  encounters  into  account.  The  Chapman  and  Cowling  conductivity 
expression  may  be  written 


<7W 


0.532  ne  e2 
(  me  kT ) % 


1 

nn  Qen 


(1) 


Several  theoretical  approaches  have  been  used  to  derive  expressions 
for  the  conductivity  of  fully  ionized  gases,  but  the  formula  derived  by 
Spitzer  and  Harm  (Ref.  20)  is  most  universally  accepted.  Because 
nn  =  0  for  fully  ionized  gases,  only  distant  encounters  are  considered. 
For  fully  ionized  gases,  Spitzer-Harm  obtain 


where 


o{  =  0.591  Ujt  e0)2 


( kT  )  \ 


in  A 


(fp  kT  )  * 
2ne 


(2) 


(3) 


is  the  ratio  of  the  Debye  length  to  an  impact  parameter  for  a  90 -deg 
deflection  of  the  electron  by  the  ion. 


The  ranges  over  which  these  conductivity  theories  are  believed  to 
provide  reasonably  accurate  results  depend  on  several  factors,  but  roughly 
the  Chapman  and  Cowling  equation  is  useful  at  temperatures  below  3000  to 
4000°K  and  the  Spitzer-Harm  equation  at  temperatures  above  10,  000  to 
12,  000°K.  No  fundamental  theory  founded  on  basic  principles  is  known  to 
exist  in  the  intermediate  range,  which  is  referred  to  hereafter  as  a  par¬ 
tially  ionized  gas.  However,  an  approximate  equation  for  the  electrical 
conductivity  of'a  plasma  having  an  arbitrary  degree  of  ionization  has  been 
proposed  by  Lin,  Resler,  and  Kantrowitz  (Ref.  8).  They  have  argued  that 
the  electrical  resistances  caused  by  neutral  particles  and  ionized  particles 
should  be  additive,  and  thus. 


This  simple  formula  appears  to  provide  reasonable  results  and  is  used 
widely  in  computations  of  a . 


Although  there  is  some  doubt  in  regard  to  defining  cross  sections  for 
charged  particles  based  on  the  concept  of  binary  collisions,  Lin,  Resler, 
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and  Kantrowitz  note  that  the  Spitzer-Harm  equation,  Eq.  (2),  assumes 
the  same  form  as  the  Chapman-Cowling  equation,  Eq.  (1),  if  one  defines 
an  effective  ion-electron  cross  section  as 


0.9  e4 

_  0n 

\  12  n 

(fo  kT)3 

(  4  77  eQ  kT  )  2 

A'  11  < 

{  e3 

2  ne 

I 

(5) 


Combination  of  Eqs.  (1)  through  (5)  results  in  an  expression  for  the  elec¬ 
trical  conductivity  for  monatomic  gases. 


a 


0.532  ne  e2 
me‘4  (kT)l/2 


(6) 


This  equation  has  been  further  generalized  to  allow  the  calculation  of 
the  conductivities  of  more  complicated  gases  and  seeded  gases.  In  a 
slightly  more  general  form,  Eq.  (6)  may  be  written. 


G 


0.532  ne  e2 
me  %  (  kT  )  l/2 


(7) 


Equation  (7),  except  for  the  numerical  coefficient,  can  also  be  derived 
by  the  so-called  "free -path  kinetic  theory"  in  which  it  is  assumed  that  the 
particles  make  instantaneous  collisions  with  each  other  but  move  freely 
between  collisions.  This  derivation  is  discussed  in  more  detail  by  Sher¬ 
man  (Ref.  21).  The  summation  term  represents  the  electron  mean  free 
path  for  gaseous  mixtures. 


Additional  terms  are  added  to  the  summation  to  account  for  collisions 
contributed  by  the  various  seed  particles  and  the  various  gas  species. 

For  example,  for  cesium-seeded  nitrogen  Eq.  (7)  would  be 


G 


0.532  ne  e2 

v  V 
me  (kT)  6 


nN 


Qn  — e  +  nCs  Qcs-e  + 


ni 


-1 


since  Qei»  given  by  Eq.  (5),  is  the  same  for  all  species.  The  numerical 
coefficients  of  Eqs.  (5)  and  (7)  were  chosen  so  that  for  the  extremes  of 
weakly  and  fully  ionized  gases,  Eq.  (7)  would  assume  the  forms  of  the 
Chapman -Cowling  and  Spitzer-Harm  equations,  respectively. 


Scalar  conductivities  (that  is,  for  B  =  0)  of  a  plasma  in  thermal 
equilibrium  have  been  calculated  from  Eq.  (7).  The  possibility  of  elec¬ 
trical  nonequilibrium  was  not  taken  into  account.  Collision  cross  sec¬ 
tions  required  in  these  calculations  and  the  methods  used  to  obtain  the 
particle  densities  are  described  in  Sections  2.  2  and  2.3. 
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2.2  COLLISION  CROSS  SECTIONS 

Much  of  the  uncertainty  associated  with  the  calculation  of  electrical 
conductivities  centers  about  the  validity  and  availability  of  collision  cross- 
sectional  data.  Brown  (Ref.  22)  and  Massey  and  Burhop  (Ref.  23)  have 
collected  much  of  the  available  cross-sectional  data.  Brown  presents 
these  data  in  the  form  of  collision  probabilities.  Also,  he  has  compared 
experimental  cross-sectional  data  obtained  by  several  investigators  in 
Fig.  1-29  of  Ref.  22  to  illustrate  the  chaotic  state  of  experimental  dis¬ 
agreement.  This  is  particularly  true  at  low  energy  levels  up  to  about 
1  ev,  which  is  the  range  of  prime  interest  for  the  a  calculation  described 
here. 

2.2.1  Electron-Ion  Cross  Sections 

By  virtue  of  the  way  that  an  effective  electron-ion  collision  cross 
section  was  found  from  the  Spitzer-Harm  equation,  Qe^  has  been  computed 
for  all  ion  species  from  Eq.  (5). 

Typical  calculated  values  of  Qe^  are  shown  in  Fig.  1  for  argon  and 
argon  seeded  with  0.  1  percent  cesium. 

2.2.2  Electron-Neutral  Cross  Sections 


Since  the  Chapman  and  Cowling  equation  (Eq.  (1))  is  derived  on  the 
basis  of  the  binary  diffusion  of  neutral  particles  and  electrons,  cross 
sections  obtained  experimentally  from  diffusion  experiments,  such  as 
the  Townsend  method,  or  calculated  diffusion  cross  sections  are  required. 

In  the  surveys  such  as  Refs.  22  and  23  it  is  not  always  clear  how  the  data 
were  obtained.  Cross  sections  obtained  from  electron  beam  or  Ramsauer- 
type  experiments  may  be  used  if  an  average  cross  section  is  computed  by 
the  method  described  in  Ref.  8  or  21. 

Experimental  electron-neutral  collision  cross  sections  for  argon 
taken  from  Refs.  22  through  26  are  shown  in  Fig.  2.  All  of  these  data 
are  believed  to  result  from  electron  diffusion  experiments,  and  the  dilem¬ 
ma  in  selecting  appropriate  values  is  apparent.  Since  the  experimental 
argon  cross  sections  vary  widely,  two  sets  of  calculations  have  been 
made  for  unseeded  argon  to  demonstrate  the  magnitude  of  possible  errors 
because  of  this  uncertainty.  The  approximate  mean  values  given  in  Ref.  22 
as  well  as  those  of  Townsend  and  Bailey  (Ref.  25)  were  used.  Lin,  Resler, 
and  Kantrowitz  (Ref.  8)  have  also  used  the  latter  values  and  have  obtained 
reasonable  agreement  with  experimental  conductivities  measured  in  a  shock 
tube. 
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In  the  computation  of  a  for  nitrogen  it  is  necessary  to  consider  both 
the  molecular  and  atomic  species.  Comparisons  of  diffusion  cross- 
sectional  data  for  molecular  nitrogen  (N2)  and  atomic  nitrogen  (N)  are 
contained  in  Fig.  3  (taken  from  Refs.  22,  23,  2  7,  and  28).  Also  shown 
are  calculated  values  of  the  diffusion  cross  section  which  were  taken 
from  Refs.  29  through  31.  The  computed  values  from  Ref.  30  result 
from  Enskog's  solution  of  the  Boltzmann  equation  and  were  obtained  by 
evaluating  the  collision  integral  for  a  Lennard- Jones  potential.  Both 
for  N2  and  N  the  analytically  determined  diffusion  cross  sections  com¬ 
pare  reasonably  well  with  experiment.  Cross  sections  selected  for  the 
present  calculations  are  those  given  by  Massey  and  Burhop  (Ref.  23), 
but  extrapolated  at  low  temperatures  as  guided  by  the  data  of  Fisk  for 
N2  and  guided  by  the  Ref.  30  calculations  for  N. 

Similarly,  experimental  cross-sectional  data  available  for  the  alkali 
metals  of  interest  as  seed  materials  are  shown  in  Fig.  4  (taken  from 
Ref.  32).  No  other  data  for  Cs,  K,  or  Na  could  be  found  in  the  literature. 
At  the  temperatures  covered  by  the  calculations  here,  it  was  necessary  to 
extrapolate  the  alkali  metal  cross-sectional  curves  to  somewhat  lower 
energy  levels. 


2.3  CALCULATION  OF  NUMBER  DENSITY 

Number  densities  of  the  various  species  have  been  obtained  by  differ¬ 
ent  methods  in  order  to  utilize  some  previous  work. 


2.3.1  Unseeded  Argon 


Since  argon  is  a  monatomic  substance  its  degree  of  ionization  at 
various  temperatures  assuming  single  ionizations  was  computed  from 


K  = 


neni 


n 


n 


/Pm  \ 
1  -  (f>2  V  kT  / 


(8) 


where  the  degree  of  ionization,  4>  (that  is,  the  fractional  number  of  ionized 
particles),  is  defined  as 


<t>  = 


ne 

nn  +  ni 


and  ne 


ni 


Equilibrium  constants  for  the  process  A  -»  A+  +  e' were  evaluated 
by  use  of  the  partition  functions  from 

2  n  me 
h2 

without  taking  into  account  the  Debye-Huckel  lowering  of  the  ionization 
potential  which  should  be  negligible  for  the  range  of  conditions  covered 
by  this  study.  Values  of  the  equilibrium  constants  for  argon  calculated 
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from  Eq.  (9),  but  multiplied  by  kT,  are  shown  in  Fig.  5.  After  the  com¬ 
putation  of  9$  at  a  given  pressure  and  temperature,  the  particle  densities 
were  found  from 


and 


ni  = 


ne 


(  0  \  Pm 

V  1  +  0  j  kT 


(10) 


nn 


/)  -  0 

V  1  +  0 


) 


Pm 

kT 


(u) 


2.3.2  Seeded  Argon 


At  temperatures  below  about  5000°K  it  would  be  sufficiently  accurate 
to  consider  only  the  ionization  of  the  seed,  but  at  higher  temperatures 
one  must  account  for  the  simultaneous  ionization  of  the  argon  atoms.  In 
the  latter  and  more  general  case,  the  derivation  of  expressions  giving  the 
particle  densities  is  straightforward  but  involves  several  algebraic  manipu¬ 
lations .  Appendix  I  contains  this  derivation,  and  only  the  results  are  given 
here. 

By  use  of  the  Saha  equation  the  degree  of  seed  ionization,  0S,  can  be 
obtained  from 


[a  _  _ 0S2  Pm  _ 

s  (1  -  0s)  (CA  +  1  +  0S)  kT  (12) 

at  temperatures  below  7000°K  when  the  degree  of  argon  ionization  is  slight. 
For  higher  temperatures  it  is  sufficiently  accurate  to  take0s  =  1.  The 
seeding  rate  factor  for  argon,  Ca*  is  related  to  seeding  rate,  S,  by 


CA  = 


(13) 


Values  of  the  equilibrium  constant,  K,  obtained  from  Eq.  (9)  for  potassium 
and  cesium  are  also  given  in  Fig.  5.  The  degree  of  ionization  of  the  argon, 
0A>  can  next  be  determined  from 


= 


0  A : 


1-0A 


Ca  Pn 


kT  [l  +  0S  +  CA  (1  +  0a) ] 


(14) 


and  then  the  various  particle  densities  expressed  in  terms  of  the  argon  and 
seed  ionization  fractions  as 


( 1  -  0a  )  CA  Pm 


kT  [  (  1  +  0S  )  +  CA  (  1  +  0A  )  ] 


(15) 
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(  ^  ~  0s  )  Pm 


kT 

[U 

+  0S  ) 

+ 

cA 

(1 

+  0 A)  3 

0 A  Ca 

Pm 

kT 

[d 

+  0s  ) 

+ 

Ca 

(l 

+  0A  )  ] 

0  s 

Pm 

kT 

[  (i 

+  0s  ) 

+ 

ca 

(1 

+  0A )  1 

(<£, 

3  +  C  A 

0. 

A  )  P 

m 

kT 

[  a 

+  0s  ) 

+ 

CA 

(1 

+  0  A  )  1 

(16) 


(17) 


(18) 


(19) 


2.3.3  Unseeded  Nitrogen 

Particle  densities  employed  in  these  calculations  were  supplied  by 
Hilsenrath  of  the  National  Bureau  of  Standards.  This  work  is  as  yet 
unpublished,  but  the  calculational  procedure  basically  follows  that  pre¬ 
viously  employed  by  Hilsenrath  in  his  air  calculations  (Ref.  33).  Some 
of  the  particle  densities  used  are  tabulated  in  Table  1. 

2.3.4  Seeded  Nitrogen 


Calculation  of  the  conductivity  of  seeded  nitrogen  also  utilized  Hil¬ 
senrath1  s  nitrogen  particle  densities.  Since  the  pressure  of  seeded  N2, 
pm,  is  equal  to  the  sum  of  the  partial  pressures  of  the  nitrogen  and  the 
seed,  the  nitrogen  particle  densities  were  evaluated  at  a  pressure  of 


Pn2  = 


cn2 

CN2  + 


1 


Pm 


The  equations  used  to  compute  the  number  densities  for  seeded  nitro¬ 
gen  are  summarized  in  Appendix  II,  but  the  general  procedure  was  to 
(1)  determine  the  N2  plasma  number  densities  at  pj^  (2)  convert  these 
to  mass  fractions;  (3)  add  the  specified  amount  of  seed  and  re-evaluate 
the  mass  fractions  of  the  new  seeded  mixture;  (4)  convert  the  seeded 
mixture  mass  fractions  to  number  densities;  and  finally,  (5)  compute  the 
degree  of  ionization  of  the  seed  by  use  of  the  Saha  equation.  The  equa¬ 
tions  are  written  for  the  addition  of  two  seed  materials  in  order  to  treat 
the  case  of  NaK.  For  computations  with  either  Cs  or  K  seed  alone,  the 
mass  of  the  second  seed  material  was  taken  as  zero. 
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3.0  CALCULATION  OF  ELECTRICAL  CONDUCTIVITY  WITH  MAGNETIC  FIELDS 
3.1  TENSOR  CONDUCTIVITY  EQUATIONS 

Since  charged  particles  move  under  the  influence  of  magnetic  fields, 
the  electrical  conductivity  of  an  ionized  gas  in  the  presence  of  a  mag¬ 
netic  field  takes  on  tensor  characteristics.  In  this  case  it  is  convenient 
to  consider  the  "effective"  directional  components  of  the  conductivity 
which  may  be  obtained  from  a  general  expression  for  the  current  density 
by  arranging  it  in  the  standard  form  of  Ohm's  law. 

With  the  following  assumptions: 

(a)  All  species  are  in  thermal  equilibrium 
(Tn  =  Tj  =  Te  =  Ts), 

(b)  Only  single  ionizations  occur. 

(c)  All  forces  are  in  equilibrium  but  the  gravitational  force 
on  the  electrons  and  electron  acceleration  is  negligible, 
and 

(d)  Electrons  are  more  free  to  spiral  than  ions  such  that 
(wr)en  >>  2(wr)in. 

Chapman  and  Cowling  (Ref.  19)  have  written  the  equations  of  motion  for 
(1)  electrons,  (2)  ions,  and  (3)  the  gas  as  a  whole.  By  combining  these 
equations  they  obtain  as  the  equation  of  charge  motion  for  a  partially 
ionized  gas  in  the  presence  of  a  magnetic  field, 

neeE'  +  (1  -  f/3)  V pe  =  Uei  +  /3/<in)BJ  +  (1  -  2f/3)  J  x  B 

+  “g— (Ken  +  [vpe  X  B  -  (JxB)  X  b]  (20) 


f  =  1  -  <t>  (23) 

and 

E'  =  E  +  V  x  B 
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In  Eq.  (20),  E'  and  B  are  the  resultant  electric  and  magnetic  field  vectors, 
respectively,  and  in  a  general  sense  may  have  components  in  all  or  any 
of  the  three  rectangular  coordinate  directions. 


Since  k  en  <<  /c^n,  <<  l,  and  “em  Brunner  (Ref.  1)  has 

simplified  Eq.  (20)  and  written  it  as 

J(1  +  b)  =  ctE  '  +  -^2 -  VPe  -  e,ere  (  J  x  -f) 


<4-4(4-i)] 


(24) 


where 


(_L_  +  -U)~l 

\  rei  ren  ) 


and 


(often  called  the  Hall  parameter) 

b  =  2f 2  rin )  (coere) 

(often  called  the  ion- slip  factor) 


In  Eq.  (24)  a  refers  to  the  scalar  conductivity,  that  is,  conductivity  with 
no  magnetic  field  present. 


By  neglecting  the  electron  pressure  gradient  and  the  last  term  on  the 
right  side  of  Eq.  (24),  Brunner  obtains  an  expressionjor  the  components 
of  current  density  parallel  to  E'  and  in  the  direction  E'  x  B  (that  is,  direc¬ 
tion  of  the  induced  Hall  field).  He  then  computed  the  effective  conductivity 
or  component  of  the  conductivity  in  the  direction  parallel  to  E'  to  illustrate 
the  Hall  and  ion-slip  effects  induced  by  a  strong  magnetic  field.  His  analy¬ 
sis  is  extended  here  to  allow  evaluation  of  all  three  main  directional  com¬ 
ponents  of  the  conductivity. 

If  the  electron  pressure  gradient  is  neglected,  Eq.  (24)  becomes 

J  (1  +  b)  =  <jE'  -  «e7e(  J  x  Hr  )  +  b  -§-  (“T  •  J  )  *25) 

The  last  Jermjwhich  was  neglected  by  Brunner  arises  from  a  term  involv¬ 
ing  (J  x  B)  x  B  which  in  a  sense  may  be  thought  of  as  the  Hall  current's. 
Hall  current.  Equation  (25)  may  be  divided  by  (1  +  b)  and  then  the  vector 
product  J  x  B  taken.  If  the  resulting  expression  for  J  x  B  is  substituted 


9 


AEDC-TDR-64-119 


back  into  Eq.  (25),  this  results  in  the  elimination  of  the  current  density 
from  the  middle  term  on  the  right-hand  side  of  Eq.  (25),  or 

J  [(1  +  b)2  +  (<yeFe)2]  =  (1  +  b)  crE '  -  <rUe7e)E'  X  -g- 


+  (were)2  +  b  (1  +  b)  ^  -|-  (-§-  ■  J  )  (26) 


Similarly,  J  can  be  eliminated  from  the  last  term  on  the  right  side 
by  taking  the  scalar  product  J  •  B  of  Eq.  (26)  and  substituting  it  back 
into  Eq.  (26).  These  operations  finally  allow  Eq.  (24)  to  be  written  in 
simple  ohmic  form  as 


(l  +  b) 


)2  +  (l  +  b): 


“1  -»•  r  oj  t  -i  /  ^ 

-  a E*  -  - = - 2_2 -  a  E’ 

L  ((yere^J  +  (1  +  b  )2  J  \ 


f  (<ye  re  )  2  +  b  (  1  +  b  )  1  I”  B  (v.  ’  .  _I_  \  I  (9.7) 

L  Kfe)!+  (1  +  b)2  J  L  B  V  B  )  J 


Equation  (2  7)  then  gives  the  components  of  the  current  density  and 
conductivity  in  the  three  directions  given  by  the  vector  terms.  The  term 
a  represents  the  scalar  conductivity  as  might  be  computed  from  the 
methods  given  in  the  previous  section,  and  the  three  components  of  the 
conductivity  are 


1. 


2. 


3. 


Parallel  to  the  resultant  E'  -field 

b 


CTeff  =r - - - 1  +  b - - 

|_(  COere  )  2  +  ( 1  +  b  )  2 


Perpendicular  to  both  E'  and  B 


„M  ,  \  - Slit - _]c 

_  (o>e  re  )  +  (1  +  b  )  2  _ 


Parallel  to  the  resultant  B-field 


P  (,COeTe  ) 2  +  b  ( l  +  b )  "1 
<reff  =  - - - - —  a 

L(WeFe)2  +  (1  +  b)  2  J 


(28) 


(2  9) 


(30) 


An  interpretation  of  these  component  equations  in  terms  of  the  effec¬ 
tive  conductivity  which  would  exist  in  various  types  of  MHD  power  gen¬ 
erators  or  accelerators  is  given  later. 
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These  results  can  be  derived  in  a  more  elegant  manner  by  writing 
Eq.  (25)  in  matrix  form  as 


Jx 

Jy  B  z  ~ 

Jz  By 

Jy 

a>ere 

+  - 

B 

J  zBx  - 

JxBz 

J  z 

JxBy  - 

JyBX 

JxBx  +  JyBy  -j-  Jz  B  z 

— 


(3 

Bx 

Ei 

By 

=  a 

E'y 

Bz 

E'z 

By  performing  the  indicated  matrix  additions  and  subtractions, 
Eq.  (31)  can  be  put  in  the  form 


An 

a  12 

ai3 

Jx 

Ei 

a2X 

a22 

a2  3 

Jy 

=  o 

Ey 

a31 

tt32 

a33 

Jz 

E  z 

and  by  inversion  of  the  matrix  |  a  | ,  Eq.  (32)  can  be  arranged  in  the 
standard  ohmic  form 

Jx  ^11  fill  fill  E  x 

Jy  =  o  fill  fin  fin  Ey  ^  ^ 

J  z  fill  fin  fin  E  z 

The  various  terms  for  both  square  matrices  are  contained  in  Appendix 
III,  Eqs.  (III-3)  and  (III-4).  With  the  matrix  terms  from  Appendix  III, 

Eq.  (27)  may  be  obtained  then  from  Eq.  (33). 

Effective  conductivities  in  practical  MHD  channels  can  be  easily 
obtained  from  Eqs.  (Ill - 3 )  and  (III-4)  in  Appendix  III.  Three  of  the  most 
z  important  cases  are  given  below  and 

illustrate  the  application  of  tensor 
■g-  conductivity.  In  each  case  one- 

^  A  y  dimensional  flow  is  considered  with 

the  orientation  shown  in  the  adjacent 

- — x  sketch.  Moreover,  induced  magnetic 

^  fields  will  be  neglected;  that  is,  low 

magnetic  Reynolds  number  is  assumed. 


3.1.1  Case  A  -  Segmented  Electrodes 

When  the  electrodes  are  segmented,  each  electrode  pair  can  assume 
a  different  base  potential  level  determined  by  the  Hall  voltage.  In  this 
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case  a  Hall  potential  exists,  but  there  is  no  Hall  current  flow.  Then, 


Bx  =  0 

Ex  =  Eh 

O 

II 

X 

By  =  B 

ii 

o 

C-H 

"C 

II 

o 

o 

II 

N 

oq 

Ez  =  Ea 

J2  =  J 

With  these  components,  Eqs.  (III-3)  and  (III-4)  reduce  to 

Jv  =  0=  (l  +  b)  Eh  +  <uere 


or 


Jy  =  0 

J  _  j  _  (1  +  b)  g  EA 

z  ~  (1  f  b)2  +  (<oere)2 


J  = 


coeTe  a  Eh 


(1  +  b)2  +  (<uere ) 2 


(1  +  b)  * 

Thus,  for  segmented  electrode  channels  the  scalar  conductivity  is  reduced 
by  the  factor  (1  +  b),  and  the  reduction  only  depends  upon  ion  slip. 


3.1.2  Case  B  —  Continuous  Electrodes 


With  infinite  continuous  electrodes,  the  Hall  potential  would  be  shorted 
out,  and  there  would  be  a  maximum  axial  Hall  current  such  that 


o 

II 

X 

DQ 

X 

II 

o 

Jx  =  Jh 

By  =  B 

Ey  =  0 

Jy  =  0 

CO 

N 

II 

o 

P3 

N 

II 

J  z  J  A 

Using  Eqs.  (III-3)  and  (III-4), 

the  current  components  in  the  x,  y,  and  z 

directions  are  found  to  be 

Jv  =  Jw  = 


((Oere  )  ea  v 


(1  +  b  )2  +  (a )ere  )2 


Jy  =  0 


J  7  J  A 


(1  +  b )  a  Ea 


(1  +  b)2  +  (a)eFe)2 

In  this  case  the  effective  electrical  conductivity  transverse  the  plasma 
is  given  by  Eq.  (28)  since  the  electric  field  is  in  that  direction,  whereas 
for  the  Hall  current  flow  the  effective  conductivity  is  given  by  Eq.  (2  9). 
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3.1.3  Case  C  -  Hall  Generator  or  Accelerator 


In  this  case  the  components 

are 

Bx  =  0 

Ex  = 

Eh 

Jx  =  Jh 

By  =  B 

Ey  = 

0 

Jy  =  0 

Bz  =  0 

Ez  = 

VB 

h  =  Ji 

and  expansion  of  the  matrix  in  Eqs.  (Ill - 3 )  and  (III - 4)  yields 

(1  +  b)  a  Ejj  c^e^e  (VB)  <7 _ 

Jx  =  Jh  =  (1  +  b)2  +  (cyer5)2  +  (1  +  b)2  +  (coe7e)2 


J 


y 


=  0 


J  J  -  (ci)ere)  g  Eh _  +  (1  +  b)  (VB)  a 

z  1  (1  +  b)2 3  +  (&je7e)2  U  +  b)2  +  (a>e7e)2 

The  effective  conductivities  for  Hall  devices  are  given  also  by  Eqs.  (28) 
and  (29). 


Equation  (2  7)  may  be  simplified  to  the  case  of  a  fully  ionized  gas  by 
taking  4>  =  1  and  hence  b  =  0;  in  which  case  the  conductivity  components 
become: 

1.  Parallel  to  the  effective  E' -field 


°eff  f 


_ I _ ]  a 

_  (ajere  )2  +  1  . 


(34) 


2.  Perpendicular  to  both  E1  and  B 


3.  Parallel  to  the  B-field 


(35) 


ffefff 


L  (coere  )2  + 


hi 1_1  a 
+  1  J 


(36) 


3.2  ION-NEUTRAL  CROSS  SECTIONS 

The  mean  free  time  between  collisions  is  related  to  the  collision  cross 
section  by 

l 

T  =  - 

v  n  Q 
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or  specifically. 


r.  =  J-  77  m‘ 
in  V  8kT 


nn  Oi 


A  /  77 

5i  ~  \ 


8kT 


1 


ne  Q  e  i 


8kT  nn  Qen 

when  the  particles  have  a  Maxwellian  distribution. 


(37) 


(38) 


(3  9) 


Determination  of  all  of  the  terms  in  Eqs.  (37)  through  (39)  except  Qin 
is  discussed  in  Sections  2.  2  and  2.  3.  However,  values  of  the  ion-neutral 
collision  cross  sections  must  now  also  be  specified. 

Some  collision  cross-sectional  data  for  cesium  and  potassium  ions  in 
argon  are  shown  in  Fig.  6a  and  for  potassium  ions  in  molecular  nitrogen 
in  Fig.  6b.  These  data  taken  from  Ref.  22  were  obtained  at  high  energy 
levels  and  require  an  extensive  extrapolation  to  the  range  of  interest  here. 
Because  the  cesium-in-argon  curve  becomes  very  steep  with  decreasing 
energy  level,  an  extrapolation  was  felt  to  be  too  impractical.  However; 
extrapolations  guided  by  other  existing  data  for  different  substances  have 
made  for  the  other  two  curves.  The  validity  of  these  gross  extrapolations 
is  certainly  subject  to  doubt.  Unfortunately,  no  existing  data  could  be 
found  for  collisions  between  potassium  ions  and  atomic  nitrogen  neutrals, 
and  consequently,  subsequent  calculations  of  the  effective  electrical  con¬ 
ductivity  in  the  primary  directions  for  fully  ionized  gases  would  be  ex¬ 
pected  to  be  less  accurate  than  those  for  partially  ionized  gases. 


4.0  NUMERICAL  COMPUTATIONS 


4.1  SCALAR  CONDUCTIVITY 


4.1.1  Argon 


Scalar  electrical  conductivities  of  unseeded  argon  at  pressures  of 
0.01,  0.  1,  1.0,  and  10  atm  were  computed  by  use  of  Eq.  (7)  by  using  the 
electron-neutral  collision  cross  sections  obtained  by  Townsend  and 
Bailey  (Ref.  25)  (shown  herein  by  the  dashed  curve  in  Fig.  2).  Additional 
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calculations  have  been  made  for  unseeded  argon  at  pressures  from  10 "2  to 
103  atm  using  the  electron-neutral  argon  cross  sections  contained  in 
Ref.  22  and  shown  in  Fig.  2  by  the  dash-dot  curve  in  order  to  make  an 
estimate  of  the  possible  errors  introduced  by  lack  of  exact  knowledge  of 
the  cross  section.  Results  of  these  conductivity  calculations  are  given 
in  Figs.  7a  and  b,  respectively,  and  are  shown  to  agree  within  about 
50  to  75  percent. 

Calculated  values  of  the  conductivity  of  argon  seeded  with  cesium  are 
given  in  Fig.  8,  whereas  Fig.  9  gives  the  results  for  potassium-seeded 
argon.  For  the  seeded  case  only  the  Townsend-Bailey  cross  sections  are 
used  to  represent  the  electron -argon  collisions.  Pressure  varying  from 
10-2  to  10  atm  is  used  as  a  parameter  for  these  conductivity  curves  to 
demonstrate  its  effect  at  seeding  rates  of  0.01,  0.  1,  1.0,  and  10  percent 
seed  weight  to  the  total  weight. 

4.1.2  Nitrogen 


Calculated  scalar  conductivities  for  unseeded  nitrogen  at  pressure 
levels  from  0.01  to  10  atm  are  shown  in  Fig.  10.  In  this  case  the  electron- 
neutral  cross  sections  used  are  those  taken  from  Ref.  23  and  shown  here 
in  Fig.  3.  In  addition,  conductivities  of  nitrogen  seeded  with  either  cesium, 
potassium,  or  NaK  in  the  amounts  of  0.  01,  0.  10,  0.  3,  0.  7,  1.  0,  and  10.  0 
percent  seed  by  weight  are  presented  in  Figs.  11  through  13,  respectively. 
In  these  figures,  seeding  rate  is  used  as  a  parameter  instead  of  pressure 
to  demonstrate  its  effect  on  the  conductivity.  Calculations  were  made  for 
NaK  consisting  of  78  percent  potassium  by  weight. 

The  initial  increase  in  the  conductivity  with  temperature  is  caused  by 
ionization  of  the  seed  material.  As  the  temperature  increases,  the  con¬ 
ductivity  reaches  an  intermediate  peak  value  where  nearly  all  of  the  seed 
is  ionized.  With  further  increase  in  temperature  the  conductivity  begins 
to  diminish.  This  occurs  because  the  curve  was  calculated  for  a  constant 
pressure;  thus  the  particle  density  and  hence  electron  density  decreases 
with  increasing  temperature  faster  than  electrons  are  created  by  N2  ioniza¬ 
tion.  The  conductivity  begins  to  increase  again  when  the  temperature  is 
high  enough  that  a  significant  number  of  electrons  are  liberated  in  the  N2« 
Conductivities  for  a  given  seeding  rate  tend  to  merge  into  those  for  the 
unseeded  case.  As  the  nitrogen  ionization  approaches  completeness, 
there  is  little  effect  of  seed  addition;  however,  at  temperatures  below 
6000°K  seeding  can  increase  the  conductivity  several  orders  of  magnitude. 
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4.2  TENSOR  CONDUCTIVITY 

Typical  values  of  the  Hall  parameter  (ueF e)  and  the  ion-slip  factor 
(b)  are  given  in  Figs.  14  and  15,  respectively,  since  their  magnitude 
indicates  the  importance  of  tensor  conductivity  effects.  They  have  been 
divided  by  B  and  b2  so  that  they  may  be  scaled  to  any  level  of  the  mag¬ 
netic  field;  the  actual  values  shown  correspond  to  B  =  1  weber/m2 
(10,  000  gauss). 

Figure  15  indicates  that  the  effective  conductivity  in  a  MHD  channel 
with  segmented  electrodes  would  be  less  than  one  percent  lower  than  the 
scalar  conductivity  at  atmospheric  pressure  and  for  seed  rates  between 
0.  1  and  1  percent  (by  weight)  when  the  magnetic  field  is  10,  000  gauss. 

The  scalar  conductivity  would  be  reduced  by  about  2  percent  if  the  B-field 
is  increased  to  20,  000  gauss.  At  lower  pressures  where  ion  slip  becomes 
much  more  important,  however,  there  will  be  a  significant  reduction  of 
the  scalar  conductivity  when  a  strong  magnetic  field  is  applied. 

To  further  illustrate  the  importance  of  tensor  conductivity  the  three 
basic  tensor  components  which  could  result  when  magnetic  fields  are 
applied  to  the  plasma  have  been  calculated  by  evaluating  aeff/  a  from 
Eqs.  (28)  through  (30).  These  calculations  have  primarily  been  made  in 
the  weakly  ionized  range  (1000°K  <  T  5  4000°K),  but  the  fully  ionized 
limit,  given  by  Eqs.  (34)  through  (36),  was  also  computed  (assuming 
T  =  15,  000°K  for  the  evaluation  of  uere  and  b).  Effects  of  the  magnetic 
field  on  the  conductivity  of  nitrogen  at  atmospheric  pressure  and  seeded 
with  0.  1  percent  potassium  are  shown  in  Fig.  16.  It  is  apparent  that  the 
effective  conductivity  in  the  direction  of  the  effective  electric  field 
decreases  rapidly  with  increasing  magnetic  field  strength.  A  similar 
comparison  is  made  in  Fig.  17  for  a  1  percent  seed  rate.  Influence  of 
seeding  rate  may  be  observed  by  comparison  of  Figs.  16  and  17  and  effects 
of  pressure  by  comparison  of  Figs.  16a  and  b  or  17a  and  b. 

The  reduction  of  scalar  conductivity  which  might  be  expected  in  chan¬ 
nels  with  continuous  electrodes  or  Hall  devices  is  indicated  by  these 
figures. 


4.3  VALIDITY  OF  CALCULATIONS 

The  accuracy  of  the  computed  conductivities  depends  of  course  upon 
the  reliability  of  the  cross  sections  used.  Various  experimental  values 
of  the  electron-neutral  cross  sections  available  in  the  literature  have 
already  been  presented  and,  at  least  for  the  case  of  unseeded  argon,  the 
magnitude  of  possible  errors  introduced  by  the  electron-neutral  term 
has  been  estimated.  In  any  event,  if  more  accurate  cross  sections 
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become  available  in  the  future,  the  present  calculations  could  be  updated. 
Additionally,  the  definition  of  the  electron-ion  cross  section  as  employed 
here  also  raises  some  questions. 

Spitzer  points  out  in  Ref.  34  that  his  theory  breaks  down  for  values 
of  Hn  A  less  than  about  five,  that  is,  at  low  temperatures  and/or  high 
pressures.  There  has  been  some  experimental  agreement  with  the 
Spitzer-Harm  theory  (Eq.  (2))  from  several  sources  (for  example  Ref.  8) 
for  values  of  £n  A  between  three  and  six.  However,  Spitzer,  himself, 
states  that  this  agreement  is  only  fortuitous.  For  many  of  the  calcula¬ 
tions  presented  here,  &n  A  falls  in  this  dubious  range.  Because  of  the 
way  the  conductivity  equation  has  been  arranged,  the  doubt  arises  in  the 
definition  of  Qei  in  Eq.  (5).  It  may  be  that  the  calculated  effective 
electron-ion  collision  cross  section  is  greater  than  the  actual  value. 

Since  the  estimated  values  of  Qei  are  from  103  to  10^  times  the  typical 
values  for  the  electron-neutral  cross  section,  the  computed  conductivity 
is  dominated  by  Qei  over  a  wider  range  than  may  actually  occur.  This  is 
particularly  true  for  seeded  gases  where  an  appreciable  electron  number 
density  exists  at  low  temperatures.  A  careful  experimental  study  of  the 
conductivity  of  seeded  gas  offers  the  possibility  of  defining  the  range  over 
which  the  Spitzer-Harm  equation  is  valid  and  still  experiment  with  gases 
at  reasonable  temperatures.  In  any  event,  one  might  question  the  calcu¬ 
lations  presented  here  when  in  A  falls  below  4.  5  to  5.  0,  and,  if  possible, 
the  calculated  results  should  be  checked  by  experiment.  Table  2  contains 
typical  values  of  in  A  for  the  present  calculations.  Calculated  conductivi¬ 
ties  of  gases  at  higher  pressure  are  certainly  doubtful;  for  example,  the 
calculation  for  argon  at  p  =  102  to  103  atm  presented  herein  should  be 
viewed  critically. 

In  the  lower  temperature  range  where  Eq.  (7)  is  dominated  by  Qen 
terms  the  calculated  results  are  believed  to  be  reliable,  since  the  Chapman- 
Cowling  equation  has  been  checked  by  several  sources  of  experimental  data. 
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APPENDIX  I 

CALCULATION  OF  PARTICLE  DENSITIES  FOR  SEEDED  ARGON 


The  mixture  pressure  of  seeded  argon  is  the  sum  of  the  partial  pres¬ 
sures  of  all  of  the  species.  If  each  specie  is  represented  by  the  perfect 
gas  law  and  the  mixture  is  in  thermal  equilibrium, 


Pm  “  ^  ni  =  (“A  +  nA  +  nA  +  n„  +  n  +  n  )  kT 


The  following  definitions  are  made: 

n  A 

*A  - 


n  *  +  n 


A. 


(I-D 


(1-2) 


<£s  = 


ns_  +  ns 


(1-3) 


A  _  c 


(1-4) 


where  n^.  and  n§  represent  the  number  densities  at  normal  conditions, 
that  is,  before  ionization  becomes  significant,  while  the  various  nj's 

(nAn>  nAp  . nSe)  are  the  number  densities  of  the  various  species 

after  ionization.  Further,  it  is  assumed  that  only  single  ionizations 
occur  so  that  nsi  =  nSe  and  =  nAe. 

Making  use  of  the  Saha  equation 


Ks  = 


2  tt  mp 


,3/2 


(kT  y 


3/2  Ze  Zf 


exp 


kT 


(1-5) 


the  equilibrium  constant  may  be  determined  from  the  right  side  of  the 
equation. 

Substitution  of  Eqs.  (I- 1)  and  (1-3)  into  (1-5)  results  in 


Ks 


cf)  p 

"s  rm 


21 


A  E  DC-TD  R-64-1 1  9 


Then  Eqs.  (1-2)  and  (1-4)  may  be  used  to  eliminate  the  remaining  number 
density  terms,  so  that 


Ks 


0  P 

_ rm _ _ 

(1  -4>s)  [1  +cf>s  +CA  (1  +^A>]  kT 


(1-6) 


In  a  similar  manner,  one  can  arrive  at  an  equation  from  which  the 
degree  of  ionization  of  argon  may  be  obtained. 


Ka 


C  .  (f>  .  p 

_ A  ^  A  _ 

(1  -<£A)  [l  +  <£s  +CA  (1  +<£A)]kT 


(1-7) 


Note  that  Eqs.  (1-6)  and  (1-7)  are  coupled.  In  principle,  the  ioniza¬ 
tion  of  the  seed  and  argon  atoms  occurs  simultaneously,  and  thus  these 
two  equations  must  be  solved  simultaneously.  Fortunately,  from  a 
practical  standpoint  it  is  not  necessary  to  do  this  since  the  seed  usually 
becomes  fully  ionized  before  the  argon  begins  to  ionize  in  any  appreciable 
amount.  For  the  range  of  the  calculations  considered  here  it  is  sufficiently 
accurate  to  neglect  the  C term  in  computations  of  from  Eq.  (1-6) 
for  temperatures  below  7000°K  since  Ca^A  <  <  Ca>  whereas  at  higher 
temperatures  the  seed  ionization  is  nearly  complete  and  a  value  of  <£s  =  1 
may  be  used. 


With  the  degree  of  ionization  for  both  the  seed  and  the  argon  deter¬ 
mined  from  Eqs.  (1-6)  and  (1-7),  the  individual  specie  number  densities 
may  be  determined  by  manipulations  of  Eq.  (1-1).  For  example,  one 
can  easily  show  that 


(1  +  <f>s)  +  Ca  (1  +  <^>A )  = 


n,  +  n  „ 
s  n  s 


+  ns 


+  n  a  +  iia  .  +  n  a 
_  A  n  A  i  A  t 


+  n  o 


Then,  using  Eq.  (1-1) 


kT  [  (  1  +  <£s  )  +  CA  (  1  +  <j&A  )  ]  (  nsn  +  ns;  ) 


1 


and  multiplying  by  nSn  gives 


( l  —  <f>  )  p 

_  _ _ _ ~  s  r  m _ 

kT  [  (  1  +  <f>s  )  +  CA  (l  +  <^>A  )  ] 

since 


(1-8) 


+  n. 


=  1  -  <f>s 


and 
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Similarly,  the  other  particle  densities  may  be  shown  to  be 


and 


_ <K  Pm _ _ 

‘Si  kT  [(l  +<f>a)  +  cA  (l  +^A)] 

CA  Pm 

An  kT  [(1  +0S)  +  CA  (1  +  0A)] 

_ C  A  Pm _ 

Ai  kT  [(1  +<f>s)  +  CA  (1  +  0A  >] 

(<£a  CA  +  Pm 

ne  =  nAi  +  ns.  =  -  [(7+<£s)  +  cA  ( 1  +  <£A  )  ] 


(1-9) 


(1-10) 


(1-11) 


(1-12) 
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APPENDIX  II 

CALCULATION  OF  PARTICLE  DENSITIES  FOR  SEEDED  NITROGEN 


Equilibrium  particle  densities  for  pure  nitrogen  which  have  been 
calculated  by  Hilsenrath  are  employed  for  all  N2  calculations.  A  few  of 
these  particle  densities  are  given  in  Table  1.  When  the  seed  is  added  to 
the  N2>  the  new  mixture  pressure  is 

Pm  =  Pg  +  Ps 

It  is  desired  to  calculate  the  conductivity  at  a  specified  value  of  pm, 
but  the  nitrogen  particle  densities  should  be  selected  at  the  N2  partial 
pressure  pg.  For  perfect  gases  in  equilibrium 

r  -  Rg  _  Pg 

^  “s  Ps 


Thus, 


p.  ■  (An)p-  <n-» 

At  pressures  given  by  Eq.  (II- 1)  and  a  specified  temperature,  the  number 
densities  or  mole  fractions  of  N2  were  found  from  tables  such  as  Table  1 
and  converted  to  specie  mass  fractions  (Yj)  by 

Yj  =  (Xj)  ( Mj )/2  Xj  Mj 


Then,  the  mass  of  seed  as  specified  by  S  is  added  to  the  mixture  and  the 
new  specie  mass  fractions  evaluated  from 


Y 


xn2  mn2  n 

2  Xj  Mj 


S) 


Y 


Nj 


XNi  MNj 
2  Xj  Mj 


(1  -  S) 
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Primed  quantities  denote  the  new  mixture.  These  in  turn  are  converted 
to  molar  fractions  or  number  fractions  for  the  new  mixture; 


Xnj  = 


(  1  -  s) 


Xn2 


Xn,  = 


2  Xj  Mj 


1  -  s 


2  Xj  Mj 


(1  -S) 


*Nj 


2  Xj  Mj 


•  1  -  s  +  _s_ 
2  Xj  Mj  Ms 


_ s_ 

M8 

1  -  s  + _ s. 

2  Xj  Mj  Ms 


Since  each  specie  is  assumed  to  be  perfect  gas,  the  total  number  of 
particles  present  may  be  determined  from 


If  two  seed  materials  are  used,  let  A  represent  the  mass  fraction  of 
one  seed  with  respect  to  the  total  amount  of  seed.  Further,  if  the  seeds 
do  not  react  with  each  other  or  with  the  N2  plasma  species,  their  original 
number  densities  may  be  similarly  computed  from 


and 


“S2  “ 


(A  )  (S  )  No  p  M 

m  m 

- KL . RT - 


(1  -  A)  S  No  Pm  Mm 

_  m  m 

Ms 2  RT 


where  the  new  mixture  molecular  weight  is  given  by 

Mm  =  (Xn2  +  Xn2.)  (28.016)  +  (Xn  +  X^)  (14.008)  +  X8't  Ms,  +  Xs'2  MSj 

Then  the  degree  of  ionization  for  each  seed  material  may  be  computed 
from 

K-.(kT)  ■  ^ -‘j.,)  °-(kT) 
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and 


Ks2UT) 


Hs2(kT) 


Finally,  from  several  of  the  above  expressions  the  number  densities  of 
each  individual  specie  may  be  evaluated  from 


n 


n 


T 


nNi  ~  ANi  nT 

nsli  =  0 Sj  nSl 

nsln  =  (  i  *“  0SX  )  nsl 
ns2.  =  0s2  r*s2 

ns2n  —  (1  ~  0s2)  ns2 

ne  “  nN2 .  +  nN i  +  nSl.  +  ns2. 

2  i  1  li  2i 
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APPENDIX  III 

TENSOR  CONDUCTIVITY  EQUATIONS 


Writing  Eq.  (25)  in  matrix  form  results  in 


Jy  B  2 

“  Jz  By 

Bx 

Ex' 

JzBx 

“  Jx  Bz 

b(jxBx  +  Jy  By  +  Jz  Bz) 

By 

E/ 

BJ 

—  o 

Jx  By 

-  J  y  B  X 

Bz 

Ez' 

(III-l) 


where  o  is  the  scalar  conductivity. 

If  the  indicated  matrix  additions  and  subtractions  are  performed  the 
resulting  equation  may  be  put  in  the  form  of 


1  +  b  -  b 


ct)e  t  e  Bz  b  Bx  By 


-(■ 


<ae  re  By 


b  Bx  By\ 


1  +  b  - 


b  B„ 


coe  re  B x  b  By  Bz 


coere  By 


-( 


+  5 


1  +  b  -  ^ 


Jx 

Jy 

Jz 


=  a 


(III-2) 


The  terms  in  the  square  matrix  are  given  as  an  to  a33  in  Eq.  (32)  of 
Section  3.1. 


Equation  (III- 2)  may  be  put  in  the  usual  form  of  Ohm's  law  by  inver¬ 
sion  of  the  first  3  by  3  matrix  term.  The  inversion  is  a  lengthy  algebraic 
process;  consequently  only  the  main  features  are  given  here.  If  Eq.  (III-2) 
is  expressed 


I A  I  IJI  =  a  lE'l 


then  by  inverting  A  it  becomes 


IJI  -  «  |  A  f  !  E  '  -  »  -iU-Li 

1  del  A 


IE1 
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The  determinant  of  A,  det  \a\,  is  found  in  the  usual  way  and  after  several  manipulations  is  simplified 
to 

det  |  A  |  =  (1  +  b)2  +  Ue7e)2 

The  adjoint  of  A,  adj  |a(,  is  evaluated  term  by  term  and  results  finally  in 


adj  |  A 


CO 

CO 


[(ct)ere)2  +  b(l  +  bl]  +  (1  + .  b) 

B*By  [(a)e7e)2  +  b(l  +  b) 

1  re  R  1 

BxB*  \u>c~cY  4.  bn  +  bll  4-  (°e  T*  R 

J 

J  B  Z 

B2  B  ” 

_A_L[(c<,e7e)s  +  b(l  +  b)]  +  ^Bz]Kl|>e7e)>  +  b(l  +  b)]  +  (1+  b)l  [MikV  +  b(l  +  b)]  BX1 

L  B  B  J  L  B2  J  J  L  B2  1  J  B 


(cje  re)2  +  b(l  4-  b)  — 


r-p-  (We)2  +  b(l  +  b)]  + 


~-\{coe~e)2  +  b(l  +  b)]  +  (1  +  b)J 


(HI- 3) 

The  terms  of  this  matrix  divided  by  the  det  |  A |  are  listed  as  through  £33  is  Eq.  (33).  In  its  final 
matrix  form  the  current  density  may  be 


Jx 

Jv 

—  ■  -  IAI 

Ex 

IT* 

y 

r  ,  ^  T -  dUJ  1  n  1 

|U  +  b)2  +  «ye7e)2J 

Jz 

E* 

z 

(III- 4) 


where  Eq.  (Ill- 3)  gives  adj  | A / .  Various  components  of  the  tensor  conductivity  in  the  most  general  case 
are  given  by  Eqs .  (Ill- 3)  and  (III-4). 
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For  the  case  of  orthogonal  fields  and  the  gas  fully  ionized,  Eq.  (III-4) 
takes  on  a  much  simpler  form.  Taking  the  fields  mutually  perpendicular 
with  the  following  directions, 

E  =  I  E 
B  =  k  B 

v  =  7  v 

Then,  Ex  =  Ez  =  0  and  Bx  =  By  =  0.  Moreover,  with  the  gas  fully  ionized 
f  =  0  and  thus  b  =  0.  Then,  Eq.  (III-4)  may  be  simplified  to 


Jx 

1 

-  c OeTe 

0 

0 

Jy 

o 

Te 

1 

0 

E ' 

1  +  (  we  re  )2 

jj 

0 

0 

1  +  ( coere f 

0 

This  is  the  same  expression  given  by  Sonnerup  (Ref.  3),  only  a 
different  orientation  was  used  here. 
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TABLE  1 


NUMBER  FRACTIONS  OF  EQUILIBRIUM  NITROGEN 


Gas 

Specie 

n2 

* 

n2+ 

T,  JK 

p  ,  atm 

p ,  atm 

10-1 

1 

10 

io-1 

1 

10 

2,000 

1.00 

1.00 

1.00 

1.80xl0"lS 

3.70xl0“19 

1.20xl0-19 

2,500 

1.00 

1.00 

1.00 

1 . 40xl0_ 1 4 

4. 30x10" 15 

1.40xl0"15 

3,000 

1.00 

1.00 

1.00 

6 .90x10” l2 

2 . 22xl0~l2 

7 . lOxlO-1 3 

3.500 

0.9992 

0.9997 

0.9999 

6.30xl0"10 

2 . 00x10“ 10 

6 . 30x10" 11 

4,000 

0 , 9944 

0.9982 

0.9994 

1.86x10"® 

5.95xl0"7 

3 . 95xl0"9 

6,000 

0.4825 

0.7925 

0.9290 

1.71x1  (T5 

1.16x10”® 

5.25xl0"6 

8.000 

0.0115 

0.0960 

0.4120 

3 . 05xl0"5 

8.33x10”® 

1.26xl0"4 

10.000 

3.90xlCT4 

4.50xl0-3 

4.60x10  2 

1.90xl0"5 

6.95x10“® 

2. 18xl0“4 

13,000 

1 . 70xl0“6 

1 . 10xl0-4 

2.00x10-3 

2.47xl0"6 

2.85x10”® 

1 . 50x10“ 4 

15,000 

1 .80x10"® 

6 . 50xl(T6 

3 . 10x10" 4 

1.90xl0"7 

8. 80x10“  6 

9 . 10x10“® 

Also  number  of  electrons  created  by  *T 


J  2. 


ionization 
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TABLE  1  (Concluded) 


Gas 

Specie 

N 

N+  * 

T.  0  K 

p  ,  atm 

p  ,  atm 

- - 

10-1 

1 

10 

10-1 

1 

10 

2,000 

2.500 

3,000 

3.500 

4 ,000 

6,000 

8,000 

10,000 

13,000 

15,000 

2 . 84xl0~ 9 
8.80x10-T 

4 . 32xlCT5 

7 . 00xl0— 4 

5.60xl0-3 

0.5180 

0.9740 

0.8543 

0.2630 

0.0540 

8 . 90x10“ 1 0 
2 . 80x10  7 

1 .38xlO"5 

2 . 30x10“ 4 

1.77X10-3 

0.2080 

0.8990 

0.9467 

0.6290 

0.3060 

2 . 85xl0— 1 0 

8 . 70xl0~ 8 

4 . 30xl(T6 
7.20x1  O’5 

5 . 60x10“ 4 

0.0710 

0.5890 

0.9375 

0.8630 

0.6740 

0.00 

0.00 

1 . 56xl0— 1 4 

1 . 15xl0-11 

1 ,88xl0“s 

1 . 03xl0— 4 

7 . 25xl0—3 

0.0727 

0.3690 

0.4730 

0.00 

0.00 

1 ,60xl0“15 

1 . 7  0x10” 12 

1 . 95x10“ 10 

1 . 70xl0— 5 

2. 18x10“ 3 

0.0243 

0. 1855 

0.3473 

0.00 

0.00 

l.eoxio-16 

1 . 18xl0“13 
1.95X10-11 

2 . 00x1 0“6 

5 . 00xl0“4 

0.0073 

0.0670 

0.1623 

Also  number  of  electrons  created  by  N  ionization. 


A  E  D  C-T  D  R-64-1 1  9 


AEDC-TDR-64-1  19 


TABLE  2 

TYPICAL  VALUES  OF  in  A  FOR  NITROGEN 


a,  S  =  0 


in  A 

T,  °  K 

p  =  0.01  atm 

p  =  1.0  atm 

p  =  10  atm 

2.000 

18.31 

17.32 

16.82 

3,000 

12.35 

11.36 

10.87 

4,000 

9.41 

6.53 

7.68 

6,000 

6.35 

5 . 56 

5.17 

8,000 

5.17 

4.20 

3.76 

10,000 

4.65 

3.56 

3.08 

13,000 

4.70 

3.15 

2.59 

15,000 

4.93 

3.12 

2.47 

•  S  -  0.01  percent  Potassium 


n  A 

T,  ° K 

p  =  0.01  atm 

p  =  1.0  atm 

p  =  10  atm 

2,000 

6.56 

5.57 

5.07 

3,000 

5.99 

4.51 

3.96 

4,000 

6.45 

4.50 

3.65 

6,000 

6.31 

5.07 

4.17 

8,000 

5.17 

4.19 

3.73 

10,000 

4.65 

3.56 

3.08 

13,000 

4.70 

3.14 

2.60 

15,000 

4.95 

3.12 

2.57 

.  S  =  1.0  percent  Potassium 


in  A 

T,  °K 

p  =  0,01  atm 

p  =  1.0  atm 

p  =  10  atm 

2 .000 

5.56 

4.58 

4.09 

3,000 

4.51 

3.46 

2.97 

4,000 

4.51 

2.98 

2.43 

6,000 

5.37 

3.27 

2.36 

8,000 

5.11 

3.76 

2.83 

10,000 

4.64 

3.50 

2.91 

13,000 

4.70 

3.14 

2.58 

15,000 

4.93 

3.12 

2.46 
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Fig.  1  Calculated  Electron-Ion  Collision  Cross  Sections  for  Argon 
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ELECTRON -NEUTRAL  COLLISION 
CROSS  SECTION.  Qen.m2 
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REF.  23 

_  _  — 

REF.  25 

— 

REF.  22 

• 

REF.  24 

o 

REF.  26 

0 


4  8  12  16  20 

TEMPERATURE,  °K 


Fig.  2  Electron-Neutral  Collision  Cross  Section  for  Argon 


2  4x10* 
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ELECTRON -NEUTRAL  COLLISION  'J,  ELECTRON-NEUTRAL  COLLISION 
CROSS  SECTION,  m2  S  CROSS  SECTION,  m2 
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uilibrium  Constant  Function  for  Selected  Elements 


ION-NEUTRAL  COLLISION  O  ION-NEUTRAL  COLLISION 

CROSS  SECTION,  m2  5  CROSS  SECTION,  m2 
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0  8  16  24  32  40  48  56  x  I03 

TEMPERATURE ,  °K 

a.  Ions  in  Argon 


x 


0  8 


16  24  32  40 

TEMPERATURE,  °K 

b.  Potassium  Ions  in  Nitrogen 


48  56  x  IQ3 

102275 


Fig.  6  Experimental  lon*Neutral  Collision  Cross  Sections  for  Various  Substances  (Ref.  22) 
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5  f  t  II  13  15 1 lO1  (  5  5  7  »  II  IS  I5»*> 

TEMPERATURE,  T ,  %  TEMPERATURE,  T ,  °K  1 - 


Fig.  9  Equilibrium  Electrical  Conductivity  of  Argon  Seeded  with 
Various  Amounts  of  Potassium 


102278 
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TEMPERATURE,  T,  *K 


TEMPERATURE,  T,  °K 


Fig.  11  Equilibrium  Electrical  Conductivity  of  Nitrogen  Seeded  with 
Various  Amounts  of  Cesium 


102280 
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electrical  conductivity. 
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Various  Amounts  of  Potassium 
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HALL  PARAMETER, 
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- Qp-O.OI  atm 


AP*0.  lOatm 
OP  -  1.00  atm 
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-  o  PARALLEL  to  E' 

-  □  PARALLEL  to  E'  x I 
“O  PARALLEL  to  B 


SYMBOLS- FULLY  IONIZED  LIMIT 


Fig.  17  Effect  of  a  Magnetic  Field  on  Electrical  Conductivity  of 
Nitrogen  Seeded  with  1.0  percent  Potassium 


